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Figure 1: Partially quenched study of the chiral behaviour of the pseudoscalar decay constant / P s [§]. 

1. Introduction 

It is a pleasure to have been invited to present lattice results on kaon physics at a conference on 
chiral dynamics. In recent years the interactions between the Lattice QCD and Chiral Perturbation 
Theory (ChPT) communities have grown very significantly, and this collaboration is very important 
for the development of our understanding of non-perturbative QCD effects in flavour physics and 
hadronic structure. On the one hand, lattice simulations are performed at values of the up and down 
quark masses which are larger than the physical ones and ChPT is used to guide the extrapolation 
of the results to the physical point. On the other hand the ability to vary the masses of the quarks 
in lattice simulations allows us to compute the low-energy constants of ChPT with unprecedented 
precision. It is to be expected that this symbiotic relation will strengthen as the precision of both 
the lattice and ChPT calculations improves further. 

The aims of this talk are i) to review recent results for important physical quantities including 
fxlfn, the form factors of semileptonic K — > % decays and the Bk parameter of K°-K° mixing; 
ii) to present a discussion of some of the theoretical conclusions, including the need to develop 
SU(2) ChPT for kaon physics and iii) to discuss future prospects including applications to K — > %% 
decays. I start however, with some brief comments about the chiral behaviour of f % . Many of the 
ideas presented below have been developed together with my colleagues from the RBC-UKQCD 
collaboration, and where appropriate I will illustrate the discussion with our results. 

2. Comments on Chiral Behaviour 

The ability to vary the masses of the quarks in lattice simulations allows us to study the chi- 
ral behaviour of physical quantities in considerable detail. Different collaborations use different 
approaches to including the NNLO (two-loop) and higher corrections, either fully or partially, cou- 
pling this with discussions of the lattice artefacts and finite volume corrections. There are a num- 
ber of talks from different collaborations at this conference discussing their procedures in some 
detail [jl|] and illustrating the dynamic activity in this field. 

In this section I will focus on one important question, whether SU(3) and/or SU(2) ChPT ad- 
equately describe the lattice data. With my colleagues from the RBC-UKQCD collaboration, we 
argue that it is SU(2) ChPT which should be used, and I now summarise the argument. Fig. |l] 
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shows the behaviour of pseudoscalar decay constant as a function of the mass of the meson. 
The right-hand upper (black) point marked with a cross corresponds to a unitary meson with 
raps ~ 420 MeV. The remaining upper (black) points have the same sea-quark mass but a vari- 
ety of valence-quark masses. Similarly in the lower (red) band the point marked with a cross is 
unitary with raps ~ 330MeV and the remaining points correspond to the same sea-quark mass but 
with different valence quark masses, keeping the meson mass below 420 MeV. The data are fit to 
NLO Partially Quenched ChPT (PQChPT) and the dashed (green) curve represents the resulting 
unitary SU(2) behaviour, which, as expected, passes through the two unitary points. The same pro- 
cedure is repeated with SU(3) PQChPT and the dot-dash (blue) curve represents the corresponding 
unitary behaviour with m s = m uc i. The important point is that the resulting value for fo (the decay 
constant in the SU(3) chiral limit) is far below the data points (60-70%), so that the validity of the 
expansion is questionable. It is for this reason that we advocate the use of SU(2) ChPT to study the 
chiral behaviour of lattice results (a further example of such a feature is shown in fig.|5] discussed in 
sec.^). Perhaps extending the SU(3) analysis beyond NLO may lead to an apparently more conver- 
gent series, and such studies are being performed but require more data at light masses and reliable 
techniques to separate chiral behaviour from small lattice subtleties. We therefore strongly prefer 
to present our results based on SU(2) ChPT and this is an area of active debate with the community. 

2.1 Kaon Chiral Perturbation Theory 

In order to apply SU(2) ChPT to kaons, the formalism must be extended. Roessl has introduced 
the corresponding Lagrangian for the interactions of kaons and pions so that he could study K% 
scattering near threshold [QJ. There are many similarities with heavy meson ChPT [Q, ||], but an 
important difference is that in the heavy quark limit ra#* = ra#, so that a B-meson can emit a soft 
pion (with coupling gBB'n) and turn into the vector meson B*. The propagation of B* mesons must 
therefore be included in diagrams. This is not the case for kaons. 

The NLO expressions for m\, fg and Bg in SU(2) ChPT (and partially quenched SU(2) ChPT) 
can be found in ref.|Q] and in section 3.2. 1 1 review the applications to decays. The applications 



of SU(2) ChPT to K —* kk decays is discussed in the talk by Hans Bijnens based on ref. [0]. 

3. The Determination of V us 

I start with a brief analysis within the standard model which I learned from my colleagues in 
the Flavianet Lattice Averaging Group (FLAG) [0]. This provides a very useful benchmark for the 
lattice results presented below. Consider the two very precise experimental results [Js[] : 

1. From the ratio of the measured rates for leptonic decays of kaons and pions we have: 

VusfK 



: 0.27599(59) , (3.1) 

Vud Jn 

where fg (f K ) is the decay constant of the K (n) and contains the QCD effects in the decay. 
Thus a precise computation of fg/ fn will give an accurate result for \Vus/Vud\- 

2. From the measurement of the differential decay rate for semileptonic K — > % decays (K^ 
decays) we obtain 

1^/40)1 =0.21661(47), (3.2) 
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Figure 2: A (preliminary) compendium of recent lattice results for /k/ f K from FLAG See [FJ] for 
publication details of the individual contributions. 

where /+(0) is one of the form factors for K — > n~£v decays at zero momentum transfer. 
Thus a precise computation of /+(0) will yield the value of \V US \. 

In addition, within the standard model we have the unitarity relation 

\V ud \ 2 + \V us \ 2 = l, (3.3) 

where \V u b\ 2 has been omitted since it is much smaller than the uncertainties on the left-hand side. 

Eqs. (|3~I|) - ([3j|) can be viewed as 3 equations for the 4 unknowns fic/fx, /+(0), V u d and V us . 
We therefore only require one additional piece of information to determine these four quantities. 
One might take, for example, the recent determination of V uc i based on 20 different superallowed 
nuclear j3 -decays ^ 

\V ud \= 0.97425(22). (3.4) 
If we accept this value of V Ll( i, then we are able to determine the remaining 3 unknowns: 

|K, | =0.22544(95), /+(0) = 0.9608(46), ^ = 1.1927(59). (3.5) 

jn 

Of course, when exploring the limits of validity of the standard model we should not assume 
the unitarity relation (fO|). I now describe lattice calculations of /k//k and /+(0). 

3.1 V us from Decays 

The difference of fy/fn from 1 is an SU(3) breaking effects, and it is important to note that 
we calculate this difference rather than the ratio itself; in the SU(3) limit we would find precisely 
1. The calculation relies on a good control of the chiral extrapolation. Recent results with two or 
three flavours of sea quarks are collected in fig.|2| The results with small quoted errors include 

^ = 1.197(3) (lf 3 ) MILC(09)g0| and ^ = 1.1 89 (2) (7) HPQCD/UKQCD (08) fllj . (3.6) 

J 71 J 71 

The preliminary FLAG summary value based on these results with Nf = 2+1 sea-quark flavours is 
Ik/ fx = 1.190(2) (10) This result is in remarkable agreement with that in eq.( p.5| ) which was 



obtained assuming the standard model and the value of V U( i in (3.4) 
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0.03 
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0.02 
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0.01152(20) 
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0.01 
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0.005 
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Table 1: The values of the form factor fo(q 
masses given in the second column Q13Q. 



2 ) 

max/ 



at the four quark masses corresponding to the four pion 



3.2 Vus from Kb Decays 

In order to determine V us from Kb decays we need a precise determination of the form fac- 
tor /n(0) = /+(0) (see eq.(3.2)). The SU(3) chiral expansion for this form factor takes the form 

/+ (0) = 1 +fi +/4 H , where f„ = 0{rrr£ K „). fi=- 0.023 contains no low energy constants and 

is well determined, fa and higher order coefficients require additional theoretical or model input 
and a benchmark value is that obtained by Leutwyler and Roos in 1984, / + (0) = 0.961 ±0.008 [|12|]. 
We see therefore that in order to be useful in extracting V us , we need to be able to compute 
/+(0) = /o(0) to better than about 1% precision. The apparently challenging target is possible 
because again it is the SU(3) breaking effects which we actually compute; the form factor is pre- 
cisely 1 in the SU(3) limit and it is the difference from 1 which we calculate. 

We start with the calculation of the form-factor fa(q% iax ), where g max = (rtiK — m n ) 2 and cor- 
responds to the kaon and pion at rest in the same frame mM. By calculating the ratio 



(K\sy 4 u\K) (K\uY4s\x) 
{n\uy A u\Ti) {K\sy A s\K) 



L/bOmax)] 



(m K + m n ) 
Am K m n 



(3.7) 



we can obtain /o(<7max) w i tri excellent precision. This is illustrated by the RBC-UKQCD re- 
sults [m]for/o(<4Jin table | 

Having obtained /o(<?max) at trie quark masses used in the simulation, we need to interpolate 
to q 2 = and extrapolate to the physical masses. Conventionally the q 2 interpolation is done by 
calculating the form factors with one of the mesons at non-zero momenta and the results presented 
below are obtained in this way. It is also possible to evaluate the form factors directly at q 2 = by 



using partially twisted boundary conditions \ \15[ |16p . 

The chiral extrapolation is more problematical and we would welcome all the guidance which 
the chiral dynamics community is able to provide. The SU(3) ChPT prediction at NLO is /o(0) = 



1 + h which both the RBC-UKQCD [JT3p and ETMC [JTVp collaborations find to lie considerably 
above their lattice results. It is not possible to exploit the two-loop results from ref. [18], since the 
results are expressed as a series with the physical pion decay constant in the expansion parameter 
and coefficients given in terms of integrals to be evaluated numerically. The NLO SU(2) ChPT 



result for /o(0) was derived in Ql9| ] and is reproduced in eq.( [3.1l| ) below, but in order to use it we 
would like results at more values of the light-quark mass to confirm that the chiral regime is indeed 
reached. The final results from the two collaborations are 



/o(0)= 0.964(5) RBC-UKQCD 0; / (0) = 0.9560(57) (62) ETMC[0]. (3.8) 
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Figure 3: Determination of V us and V ul i based on the lattice results for Vus/Vud an d the RBC-UKQCD 
result for V us . The vertical blue band corresponds to V uc i in [gp. The black curve is the unitarity condition 
|V„ d | 2 + |V,v| 2 = l. 



These values are lower than those obtained using analytical methods based on the two-loop ChPT 



formulae of [18]. 



Comparing the lattice results in eqs.(3.6) and (3.8) with the expectations from the standard 
model ( |3.5[ ) (which, it should be noted, was based on the value of Vm in (3.4)) we see that there is 
very little room for new physics contributions to the violation of the unitarity of the first row of the 
CKM matrix. This is illustrated in fig.|3] 1 and will be quantified in [|7j|. 

3.2.1 SU(2) Chiral Perturbation Theory for K& Decays 

In this section I review the applications of SU(2) chiral perturbation theory to K^ decays at 
q 2 = <7^ ax where it is natural and at q 2 = where, because the pion is hard, it is not [jlSfl. 

Kn decays at ^ax- SU(2) chiral perturbation theory can naturally be applied to Kb decays near 
the end point of phase space where the momentum of the pion is small (in the rest frame of the 
kaon). Although we cannot perform simulations at the unphysical Callan-Treiman point where 
q 2 = m 2 K — m 2 n , as m u = mj — > the value of q 2 at the edge of the physical phase space, q 2 ^ = 
(mg — m K ) 2 approaches the Callan-Treiman point and we have the relation 

m 

/oteLO — n V- 1 - 26 ' (3 - 9) 

where p K > and / are the kaon and pion decays constants in the SU(2) chiral limit, and the value 1.26 
is estimated on the basis of the chiral behaviour of the ratio of decay constants in the RBC-UKQCD 
simulation [§] 2 . On the other hand when we look in table| at the RBC-UKQCD values of /o^Lx) 



obtained with pion masses ranging from 670 MeV to 330 MeV J13Q there is little indication of 
convergence towards 1.26 in the chiral limit. We investigated this using SU (2) ChPT and found that 
although the chiral logarithms are of approximately the correct size to account for the difference 



1 1 am grateful to Andreas Juttner for providing this figure and to the Flavianet Kaon Working Group for developing 
this way of presenting the results. 

2 Thanks to Enno Scholz for providing this number. 
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between the results in table [T] and 1.26 they have the wrong sign, tending to make the form factor 
decrease towards the chiral limit. The result is 



f ( 2 ^ f { 



1 max ) 



f 



11 mi ml X\ %2 ? 

—r y. — 777 l°g TT + 1 — 7 m x + 71 — 77? m % + ' ' 

4 (4tt/) 2 & A| 4tt/ (4tt/) 2 



(3.10) 



The coefficients Ai 2 cannot be evaluated within SU (2) ChPT They can be estimated however, by 
converting the results from SU (3) chiral perturbation theory, where there are no unknown constants 
at one-loop level except for f( K > / f [ j20| , 21]. This conversion suggests that they have the correct 



sign and magnitude to account for the difference of the observed results from 1 .26. 

It should be noted that this rapid chiral behaviour at the edge of phase space is not limited to 
the form factor of the kaon, but also applies to the heavy mesons where 

JO (<7max) and Jo (<7max) ► ~ T" j 

m K ->0 J m x ^0 J 

and yet in the region of quark masses where we have data, the /o(g max ) are below the expected 
ratios of decay constants. As an example consider fig. 14 in ref.[[22|] where the form factor fo for 
B — > k semileptonic decays does not appear to approach the expected value of // as q 2 — > g max . 

Kb decays at q 2 = 0. At first sight it appears that SU(2) ChPT should not be applied at q 2 = 0, 
since in this case the final-state pion is hard; in the rest frame of the kaon its energy is m^/2. 
However, we argue in [|l^] that the chiral logarithms come from soft internal loops and that it is 
possible to evaluate them: 

A(o)=/ + (o) = F +( 1 -ii^i°s@)+< + »4) o.ii) 

^^-Hw'^H"' 2 *)- <3i2) 

where F + _ and c + _ are unknown constants which depend on m s but not on m u d. In arriving at 



( J3- 1 1[ ) and ( |3.12| ) we use integration by parts and equations of motion to reduce the infinite chain of 
higher order operators which contain additional derivatives on the hard external pion to the leading 
operator. It may be possible to formulate this approach to ChPT with hard external pions in terms 
of an effective theory in which hard and soft pions are separated, but this has not yet been achieved. 

The approach described here can also be applied to the semileptonic form factors at other 
values of q 2 and to other processes. It has been applied to nonleptonic kaon decays [Q] and in the 
preceding talk Bijnens has referred to the approach as Heavy Pion Chiral Perturbation Theory. 

Since the chiral extrapolation is the major source of systematic uncertainty for the lattice de- 
termination of V us from K(3 decays, it is important to have all the possible theoretical information 



to guide us, and in particular it would be very helpful to extend the result in ( p.ll| ) to NNLO. 
4. B K 

Calculations of the parameter Bk which contains the non-perturbative QCD effects in K° -K° 
mixing have been performed for about 20 years. Bk is defined in terms of the matrix element of 
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Figure 4: Recent unquenched results for Bk together with the average proposed by P.Boyle at the 2009 
Kaon International Conference. 



the single AS = 2 operator which contributes in the Standard Model: 

(K° | (jy M (l - y 5 )d) {sf{\ - n)d) \K°) = ^f K m 2 K B K (^) . (4.1) 



If the lattice formulation has good chiral and flavour symmetries then the AS=2 operator in (11) 
renormalizes multiplicatively. The renormalization can be performed nonperturbatively. Recent 
unquenched results for the scale invariant parameter Bk are shown in fig.^. I highlight two 
calculations with results which are in remarkable agreement. The first is the RBC-UKQCD re- 
sult presented in [g g], B K = 0.720(13)(37) which was obtained from its domain wall fermion 
(DWF) datasets at a single lattice spacing but with an estimate of 4% for the discretisation uncer- 
tainties; this is the largest component of the quoted error. (Bk = 0.720(13)(37) corresponds to 
B^ s (2GeV) = 0.524(10)(28).) The second calculation by Aubin, Laiho and Van de Water ob- 
tained Bk = 0.724(8) (28) [ p4| ] and was obtained using domain wall valence quarks and staggered 
sea quarks at two lattice spacings. The taste unitarity violations were removed using SU(3) CHPT 
at NLO (including some of the NNLO terms). In both cases the AS = 2 operator was renormalized 
non-perturbatively using the same techniques. 
Given these results I am happy to take 

B K = 0.722(40) (4.2) 

as the current best value as was proposed by RBoyle at the 2009 Kaon International Conference. 



5. K^nn Decays 

An understanding of the longstanding puzzle of the AI = 1/2 rule and the theoretical deter- 
mination of e'/e, whose experimental measurement confirmed the existence of direct CP- violation 
in kaon decays, are major long-term goals for the lattice community. In 2001, two collaborations 
published some very interesting quenched results on non-leptonic kaon decays, in particular, 



Collaboration(s) 


Re A /Re A 2 
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Experiments 


25.3 ±1.8 

9-rl2 

22.2 


-(4.0 ±2.3) x 10~ 4 

(-7h--2)x10~ 4 
(17.2±1.8) x 10~ 4 
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Figure 5: The mass dependence of the K — > % matrix element. m\, m x and m z are the masses of the sea and 
valence light quark and the valence "heavy" (strange) quarks respectively. 

where Aq and A2 are the decays into isospin 1 = and 1 = 2 two-pion final states. These results 
were obtained not only in the quenched approximation, but also at lowest order in the SU (3) chiral 
expansion and at relatively heavy masses {m n > 400 — 500 MeV). At Lowest Order in SU(3) ChPT 
the K — > 7171 decay amplitude is given in terms of K — > 71 and K — * vacuum matrix elements, i.e. with 
at most one hadron in the initial and final state. In spite of the limitations of these calculations, the 
authors did achieve the control of the ultraviolet problem, i.e. the numerical subtraction of power 
divergences and the renormalization of the weak operators. This is highly non-trivial. 

The RBC/UKQCD collaboration have repeated the calculation with the 24 3 DWF ensembles 



in the pion-mass range 240-420 MeV Q25 , [260 ■ F° r illustration consider the determination of a 2 7> 
the LO LEC for the AI = 3/2 (27,1) operator: 



0(27,1) = {sd) L {(uu) L -(dd) L } + {su) L (ud) L , 



(5.1) 



where the subscript L stands for left. Satisfactory fits for the mass dependence were obtained using 
NLO SU(3) ChPT, but again the corrections were found to be very large, casting serious doubt 
on the approach. This is illustrated in fig.|5[ the dashed curve represents the curve in the SU(3) 
limit (m u = = m s ) and the value of this curve in the chiral limit is much below the data points, 
demonstrating that the one-loop corrections are very large. Thus the use of soft-pion theorems is 
not sufficiently reliable and K — > %% matrix elements have to be computed directly. 

In general the computation of matrix elements with multi-hadron states is not understood theo- 
retically. For K — * 7171 decays however, under the assumption that the s-wave phase-shift dominates 
the rescattering and that inelastic states can be neglected the theoretical framework is understood. 
Consider first the propagation of two pions, between time and t. Four diagrams contribute to the 
two-pion propagator: 




t 
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Figure 6: The correlation function and effective mass for the (itn)i=2 propagator (top line) and for the 
{nn)i = Q propagator (middle line) as a function of time t. The bottom line shows the (tztz)j = q propagator with 
the contribution from diagram V removed. The data is from an exploratory calculation by RBC-UKQCD 
(courtesy of Qi Liu). 

If the two pions are in an 1=2 state then only diagrams D and C contribute and the correlation 
function is proportional to D-C. For the propagation of the 1=0 state all four diagrams contribute 
and the correlator is proportional to 2D+C-6R+3V. The principle difficulty in the evaluation of 
K — > (7F7f)/ = o decay amplitudes is the evaluation of diagrams with a vacuum contribution such as 
diagram V in the propagation of two pions. To illustrate this point consider the plots shown in 
fig. |6| which show the correlation functions for two pions (left-hand panels) and the corresponding 
effective masses (right-hand panels), both as functions of the time t. The data is from an exploratory 
study by the RBC-UKQCD collaboration and I am grateful to Qi Liu for providing them. The top 
row contains the results for the 1 = 2 state and as expected these are very clean. The middle row 
on the other hand, shows the corresponding results for the 1 = state and we see that large errors 
start at an early time. The bottom row shows the 1 = correlation function but with the V diagram 
removed, highlighting that the problem is largely due to the precision with which we can do the 
vacuum subtraction from the correlation function. In order to make further progress we need to try 
to develop techniques which would allow us to extend the range of t where the errors are small and 
to use this range most effectively. Attempts to do this are under way. 

6. Conclusions and Prospects 

These are exciting times. Unquenched lattice simulations are presenting results with good 
precision with very light quarks and much effort is now being devoted to considering how lattice 
calculations will best help to clarify any new physics discovered at the LHC. The consistency of 
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results with different actions is impressive and important and adds hugely to our confidence. The 
chiral regime is being mapped out for the spectrum and decay constants with some simulations 
bring performed with light quarks close to the physical ones. The more expensive simulations 
with good chiral and flavour properties, such as the domain wall fermions which were used by 
RBC-UKQCD in obtaining many of the results presented above, will be used to extend the range 
of physical quantities which can be studied in lattice kaon physics. I imagine that by the next 
Chiral Dynamics Workshop in 2012, the chiral behaviour of many fundamental quantities will be 
well understood and the low energy constants will be determined with excellent precision. For this 
to be achieved, there will have to be continued close collaboration between the ChPT and Lattice 
communities. 

In this talk I have not had the opportunity to discuss how we perform the renormalization 
of bare lattice quantities or operators. I want to stress however, that the precision of lattice cal- 
culations has now reached the point where we need significant cooperation from the perturbative 
QCD community. Many (but not all) lattice calculations are implemented using non-perturbative 
renormalization, in which the renormalization conditions are performed non-perturbatively and 
perturbation theory using the lattice action is completely avoided. For this to be possible the renor- 
malization conditions have to be ones which can be simulated, such as ones based on the evaluation 
of Green functions at external momenta which serve as the renormalization scale [|7|] or on the use 



of the Schrodinger functional [28]. Higher order perturbative calculations on the other hand tend 
to be performed using dimensional regularization which cannot be simulated. Thus, in order to 
combine lattice results for operator matrix elements with perturbatively computed Wilson coeffi- 
cient functions, a matching between two continuum schemes must be performed in perturbation 
theory. The matching factor is frequently only available at one-loop order, leading to a significant 
error by today's standard (an important example is the evaluation of m s with domain wall fermions 



as discussed in Q29Q, where the error in the determined mass in the MS scheme is largely due to 
this continuum calculation.) Cooperation between the perturbative QCD and lattice communities 
to define renormalization schemes which can be implemented in lattice calculations and yet which 
are convenient for higher order perturbative calculations is now necessary. 

I hope that in this talk I have demonstrated that a huge amount has already been achieved in 
lattice kaon physics, but that major challenges still lie ahead. 

Acknowledgements I warmly thank my colleagues from the RBC-UKQCD and FLAG Collab- 
orations with whom I have developed much of my understanding of the material of the talk. I 
acknowledge partial support from STFC Grant ST/G000557/1 and EU contract MRTN-CT-2006- 
035482 (Flavianet). 
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